Spontaneous Formation of Dipolar Metal Nanoclusters 
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The adsorption of Ags and Ag4 clusters on the a-Al2O3(0001) surface is explored with den- 
sity functional theory. Within each adsorbed cluster, two different cluster-surface interactions are 
present. We find that silver clusters simultaneously form both ionic bonds with surface oxygen 
and intermetallic bonds with surface aluminum. The simultaneous formation of disparate electronic 
structure motifs within a single metal nanoparticle is termed a "dipolar nanocluster" . This coexis- 
tence is ascribed to the similar bond enthalpies of Ag-Al and Ag-0 bonds, and its importance for 
nanoparticle catalysis is highlighted. 



There is great fundamental interest in understanding 
how transition metals and oxides are affected by con- 
tact with each other. Since transition metal/oxide inter- 
faces exert significant controllable influence on material 
properties, many current and potential applications rely 
on these heterostructures, including catalysts for auto- 
motive pollution control [l| and fuel cells 0], as well 
as nanoscale biosensors. Jp Recent theoretical ^ |E IS 
and experimental 0, H, Q studies of size-selected nan- 
ocluster deposition onto oxides highlight size depen- 
dence of the chemical and physical properties. Fur- 
thermore, numerous experimental and theoretical inves- 
tigations emphasize the role of the oxide support in 
changing the catalytic ability of these metal-oxide sys- 
tems. [IE 111 El El 111 III 111 113 

The combination of Ag and AI2O3 is special, because 
of the close competition between ionic (Ag-0) and in- 
termetallic (Ag-Al) bonding. Ag and Au have similarly 
low oxide formation energies, [l^ However, Au bonds to 
Al much more strongly than to O, whereas the bonds 
that Ag makes to Al and to O (Ag-Alwl.gi eV and 
Ag-0~2.28 eV) are more similar than any other element 
(only Cu is close). ,19] 

In the present Letter, we report a novel consequence 
of this bonding competition: small supported Ag„ clus- 
ters exhibit two coexisting structural and electronic re- 
lationships to an a-Al2O3(0001) (a-alumina) substrate. 
The proximity of such different states within one clus- 
ter is fundamentally interesting, and has ramifications 
for understanding and improving noble metal nanoclus- 
ter catalysis. We use first-principles density functional 
theory (DFT) to study the bonding of three- and four- 
atom Ag clusters to the Al-terminated (0001) surface of 
Q;-Al203. We find that bonding competition strongly in- 
fluences the stable cluster adsorption geometries, and is 
directly responsible for inducing unusual electronic states 
in these clusters. 

DFT calculations were performed with a generalized- 
gradient approximation exchange-correlation func- 
tional. Ill] Geometry optimizations were carried out 
using an in-house code, and calculation of Born effective 
charges 21] and orbital-projected density of states 
(PDOS) were done using the ABINIT software pack- 
age. [23] All calculations were converged using a 2 x 2 x 1 



grid of Monkhorst-Pack /c-points. [23] Norm-conserving 
optimized pseudopotentials l24l with the designed 
nonlocal method for metals |25l |2^ were constructed 
using the OPIUM pseudopotential package. The 
Kohn-Sham orbitals are expanded in a plane-wave basis 
set truncated at 50 Ry. 

The a-Al203 surface is modeled by a slab geometry 
supercell with an in- plane (\/3 x V3)m0° unit ceU and 
periodic boundary conditions. The slabs consist of flve 
AI3O9AI3 tri-layers, making the surfaces Al-terminated. 
At least 12 A of vacuum separate periodic images in the 
(0001) direction. The theoretical AI2O3 in-plane lattice 
constant of 4.798 A was used (4.759 A experimental '53]). 

At equilibrium, the surface Al layer is only z « 0.1 A 
above the O layer. In-plane relaxation yields nearest- 
neighbor 0-0 distances from 2.63 A to 2.94 A. Structural 
details of the alumina surface are consistent with other 
modeling studies p^ . 

Chemisorption of Ag on the g-AloOs (0001) surface 
strongly favors cluster formation |30l Isil ] because of 4- 
5% mismatch between the surface 0-0 distance (2.76 A) 
and the Ag-Ag distance (2.88 A) in bulk Ag. To find the 
minimum-energy structures for Ags, for planar Ag4, and 
for pyramidal Ag4, each cluster was initially placed three 
different ways. Interfacial Ag atoms were started at top, 
bridge, or hollow sites of the oxygen lattice, and each 
system was allowed to relax fully. In each calculation, 
the top two surface tri-layers of alumina were allowed to 
relax in all directions, the third tri-layer was relaxed in 
the direction normal to the surface plane, and the bot- 
tom two tri-layers were constrained to their bulk alumina 
coordinates. 

For all starting positions, geometry optimization of 
Ag3/Al203 leads to the same final structure (Figure^, 
with all Ag atoms in hollow sites of the top O layer lat- 
tice, one above a surface Al atom, and two above subsur- 
face Al atoms. The cluster chemisorption energy i?ads is 
2.09 eV. The cluster tilts by 30° with respect to the sur- 
face plane, with the Ag above the surface Al farthest from 
the surface. Defining z as the distance of Al atoms from 
the topmost oxygen layer in the surface normal direction, 
the values in bare AI2O3 are 0.096 A for the topmost Al 
of the unit cell. In response to cluster adsorption, the 
nearest surface Al atoms move vertically by Az =-0.47, - 



2 




(a) 



(b) 



FIG. 1: Induced charge density Ap diagrams for the optimized 
structure of Ags/A^Oa. O, Ag (connected), and Al are shown 
as black spheres of decreasing size. Adsorption causes electron 
flow from dark to light regions. Topmost Ag and Al atoms are 
labeled. Iso-surface values are ±0.02 e^/A^ (a) Side view, 
(b) Top view. 



0.55, and +0.57 A. One Ag and one Al atom shift upward 
together. The inward relaxation of the the other nearby 
surface Al atoms induced by metal cluster adsorption has 
been observed in other theoretical studies [s^ . 

The relaxed interatomic distances strongly suggest two 
different bonding motifs for the Ags cluster. The raised 
Ag and Al make a short 2.65 A bond, which is just about 
the sum of their covalent radii (rAg=1.53 A, rAi=1.18 A). 
The shortest distance from this Ag to O is 3.66 A, much 
longer than the sum of their ionic radii (rAg+=1.14 A, 
ro2-=1.24 A). The other two Ag atoms have short 
«2.46 A distances to O, quite in line with ionic bonding. 
These Ag atoms also have long distances of 2.95 A to sub- 
surface Al, suggesting little if any covalent interaction in 
that case. So the bond length data can be summarized 
as a covalent (intermetallic, IM) bond between the raised 
Ag and Al, and ionic Ag-0 bonds for the other two Ag 
atoms. 

We study Ap, the change in charge density induced by 
the cluster adsorption, to visualize the electronic cluster- 
surface interactions. The side and top views of the iso- 
surfaces for Ap are shown in Figure ^ with electronic 
charge flowing from dark to light regions. Figure ^ shows 
gain of charge between the raised Ag and Al, indicative 
of Ag-Al IM bond formation. The other two Ag atoms 
show a significant loss of electrons, with a corresponding 
gain of electrons for the nearest topmost surface oxygens. 
Therefore, these interactions are chiefly ionic bonds. 

This electronic description of Ags adsorption also pro- 
vides insight into the observed surface relaxation. The 
ionic bonds formed between Ag and O leave these O 
atoms less capable of bond ing to surface Al atoms, using 
a bond-valence argument. |33 Therefore, these adjacent 
surface Al atoms relax inward (below the top O layer). 



TABLE L Principal values of Born effective charge tensors 
and electronic rf-band centers e^j for supported cluster Ag 
atoms. For each atom, the largest principal value of each 
sign, along with the corresponding principal direction, are re- 
ported. The principal directions (shown in Figure ^ in all 
cases point inward toward the cluster center and along the 
surface normal. The angle 6 between the principal direction 
and the surface normal are tabulated. 





Principal Value 


e; 


Ed, eV 


Ag3 Ionic 


1.11 


32 


-3.42 


Ag3 IM 


-0.23 


60 


-4.32 


Boat ionic 


1.45 


31 


-3.14 


Boat IM 


-0.42 


41 


-3.66 


Boat IM 


0.11 


87 


-3.66 


Candlestick ionic 


1.22 


26 


-3.63 


Candlestick IM 


-0.35 


59 


-4.56 


Pyr ionic 


1.48 


23 


-3.39 


Pyr IM 


-0.27 





-4.16 


Pyr IM 


0.36 


85 


-4.16 


Top Ag 


-0.11 


isotropic 


-3.78 



so they can form bonds with subsurface oxygen. The 
increase in electron density near the IM Al reduces its 
electrostatic interaction with oxygen and results in the 
observed outward relaxation. 

We propose referring to these clusters as "(electric) 
dipolar nanoparticles." (Magnetic) dipolar nanoparticles 
have been reported |33|, but we know of no previous re- 
port showing spontaneous formation of electric dipoles 
on nano-sized supported metal particles. 

To make precise the magnitude and orientation of the 
nanoparticle dipole, we apply the modern theory of po- 
larization. [21I The Born effective charge tensor Z* is 
found [3^ l36| . where the tensor element Z*^j gives the 
change in electric polarization component Pi as atom 
a moves along direction j. There is no requirement 
that this mixed second derivative tensor be symmet- 
ric, and the tensors of the supported clusters are highly 
anisotropic. We have obtained the principal values of 
the charge tensors and for each supported Ag atom the 
largest principal values (of each sign, where applicable) 
and their corresponding principal directions are reported 
in Table n 

The dynamical charge values support the above inter- 
pretation of the Ag/Al203 interactions: ionically bound 
Ag atoms lose charge through interaction with surface O, 
resulting in positive principal values. The IM Ag atom 
gains charge from the bonding with surface Al, result- 
ing in a modest negative Z*, but most of the charge is 
shared, not closely associated with Ag motion. In some 
cases (boat, pyr) the IM atoms have a more complicated 
Born effective charge tensor, with positive and negative 
principal values. The principal directions all point in- 
ward to the center of the cluster and out of the surface 
plane (Figure OJ. Increased negative charge on the IM Al 
atom and ionic O atoms in the adsorbed geometry rela- 
tive to the bare surface compensate for the net positive 
charge localized on the cluster. 
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FIG. 2: Density of states projected onto the s-orbital of Ag 
atoms in Aga/AbOa- 



Ionic and IM Ag-surface interactions are also revealed 
in PDOS analysis. Covalent IM bonding causes orbitals 
to mix, leading to intensity in the IM Ag s-PDOS below 
the Fermi level (Figure ^ . The s orbital of the ionic Ag 
shows less intensity at bonding levels, and is dominated 
by substantial intensity above the Fermi level. The frac- 
tional fillings of the IM and ionic Ag s orbitals are 0.69 
and 0.35, respectively. 

The projection onto Ag d-statcs furhtcr demonstrates 
how the coexistence of two cluster/surface interactions 
affects the reactivity of the supported Ag atoms. Our 
analysis shows that the filling of the supported Ag ri- 
bands is constant and near unity for all Ag atoms. Tabled 
lists the average energy of the rf-band projections (cd) of 
supported Ag atoms with respect to the Fermi level, 
is well established as a predictive parameter for assessing 
reactivity 37]. The 0.90 eV shift in between the ionic 
and IM Ag atoms in Ag3/Al203 is larger then what can 
be achieved through perturbations such as strain jsl] and 
is more in line with the extent of shift brought about by 
significant coordination change j39| or introduction of a 
metal hetero metal atom to a surface |40j . 

Having built a model of the chemisorption bonding in 
the Ag3/Al203 system, we examine how the chemisorp- 
tion of Ag4 onto AI2O3 expands these ideas. As with the 
Ag3/Al203 system, calculations of Z* and PDOS were 
carried out on all optimized cluster geometries, and key 
results are discussed. 

Pyramidal Ag4 ("pyr") bonds to AI2O3 
(-£'ads=l-96 eV) such that the three base atoms 
tilt and interact with the surface very similarly to 
Ag3/Al203. This and near-zero principal values of Z* 
for the pyramidal top Ag atom imply that the base 
screens the top atom from electrostatic interaction with 
the surface, consistent with previous results concerning 
the length scale of metal-oxide interactions. A single 
Ag-Al IM bond is formed, and the other two Ag atoms 
participate in ionic interactions with the topmost surface 
oxygen. The Al Az values are similar to those found 
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FIG. 3: Top views of optimized cluster geometries. O, Ag 
(connected), and Al are shown as black spheres of decreasing 
size. Unit vectors of the principal directions of the Born ef- 
fective charge principal values are indicated by arrows. Only 
the topmost nine O and three Al surface atoms of the unit 
cell are shown, (a) Aga (b) Pyr (c) Boat (d) Candlestick 



in the Ag3/Al203 structure. The top Ag atom of the 
pyramid in the supported geometry is not in the Ag 
trimer hollow, but is shifted almost to the bridge site, 
above the two charge-depleted ionic Ag atoms. 

The bonding competition between Ag-Al and Ag-0 
directly leads to two metastable minima for the planar 
Ag4/Al203 system. Each starting structure is a parallel- 
ogram parallel to the surface. The cluster buckles signif- 
icantly as it chemisorbs. We find two local minima, each 
with all Ag atoms in hollow sites of the O lattice; the 
structures differ in their registry with the top Al sublat- 
tice. The ground state ("boat") has two Ag atoms above 
surface Al, while the other local minimum ( "candlestick" ) 
has only one Ag-Al interaction. 

The boat and candlestick Ag4 clusters have E'ads values 
of 2.14 eV and 1.49 eV, respectively. The boat is ener- 
getically favored, suggesting that a balance of IM and 
ionic bonding stabilized both. In fact, the boat exhibits 
larger positive and negative Z* values than the candle- 
stick. Therefore, the energetically favored Ag4 boat clus- 
ter will exhibit even stronger dipolar nanoparticle prop- 
erties than Ag3. 

In conclusion, we find that bonding competition be- 
tween Ag-Al and Ag-O gives rise to ionic Ag-0 and in- 
termetallic Ag-Al interactions between Ag cluster atoms 
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and the alumina surface. The proximal coexistence of 
these interactions results in the formation of dipolar 
nanoparticles. The electronic and structural effects are 
closely related, with IM and ionic Ag-surface bonding 
favoring outward and inward Al motion, respectively. 
We find consistent results and interpretations of induced 
charge density, Born effective charges, and projected den- 
sity of states in all four optimized cluster geometries 
(Ag3, pyr, boat, and candlestick). Principal values of Z* 
show that ionic Ag atoms with positive charge and IM 
atoms with negative charge can be clearly distinguished, 
while stacked Ag atoms are mostly screened. The coexis- 



tence of ionic and IM bonding in supported clusters may 
be possible to create in other metal-oxide combinations, 
where the oxide and IM bond enthalpies are inherently 
similar or by tuning the competition by means of surface 
modification. 
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